Bimodality - a Sign of Critical Behavior in Nuclear Reactions 
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The recently discovered coexistence of multifragmentation and residue production for the same 
total transverse energy of light charged particles, which has been dubbed bimodality like it has been 
introduced in the framework of equilibrium thermodynamics, can be well reproduced in numerical 
simulations of the heavy ion reactions. A detailed analysis shows that fluctuations (introduced 
by elementary nucleon-nucleon collisions) determine which of the exit states is realized. Thus, we 
can identify bifurcation in heavy ion reactions as a critical phenomenon. Also the scaling of the 
coexistence region with beam energy is well reproduced in these results from the QMD simulation 
program. 

PACS numbers: 24.10.Lx, 24.60.Lz, 25.70.Pq 



Recently, the INDRA collaboration has discovered [l| 
that, in collisions of heavy ions - Xe+Sn and Au+Au be- 
tween 60 and \QQA.MeV incident energy -, in a small in- 
terval of the total transverse energy of light charged par- 
ticles {Z < 2), E± 12 , a quantity which is usually consid- 
ered as a good measure for the centrality of the reaction, 
two distinct reaction scenarios exist. In this E± 12 inter- 
val, in forward direction - i.e. quasi-projectile ~, either a 
heavy residue is formed which emits light charged parti- 
cles only, or the system fragments into several intermedi- 
ate mass fragments. This phenomenon has been named 
"bimodality". In addition, as shown in the mean 
E±_ 12 value of this transition interval scales with the 
projectile energy in the center of mass of the system for 
Au-f Au reactions, between QOA.MeV and 150 A.MeV. 

This observation has created a lot of attention, be- 
cause a couple of years before, the theory has predicted 
0, [1] that in finite size systems, whose infinite counter- 
parts show a first order transition, the system can - for 
a given temperature - be in either of the two phases if 
this temperature is close to that of the phase transition. 
Assuming that E± 12 is a measure for the excitation en- 
ergy, and acts as the control parameter of the system, it 
is tempting to identify the residue with the liquid phase 
of nuclear matter, and the creation of several medium or 
small size fragments with the gas phase. The experimen- 
tal observation would then just be a realization of the 
theoretical prediction. 

If this were the case, the longstanding problem to iden- 
tify the reaction mechanism which leads to multifrag- 
mentation would be solved. This problem arrived be- 
cause many obscrvablcs could be equally well described 
in thermodynaniical or statistical theories 0, Q a-s in dy- 
namical models 0, 01 , although the underlying reaction 
mechanism was quite different. The statistical models 
assume that at freeze out, when the system is well below 
normal nuclear matter density, the fragment distribution 
is determined by phase space. 



In dynamical models, on the contrary, fragments are 
surviving initial state correlations which have not been 
destroyed during the reaction, and equilibrium is not es- 
tablished during the reaction. A detailed discussion of 
how the reaction proceeds in these models can be found 
ini. 

To quantify the bimodality, one may define as in ref. m 

0-2 {Zmax ^ Z„iax-l) / {Zmax + Z^ax-l) (1) 

where Zmax is the charge of the largest fragment, while 
Zmax-i is the charge of the second largest fragment, both 
observed in the same event in the forward hemisphere - 
at polar angles 9cm < 90° - in the center of mass of the 
system. If the system shows bimodality, we will observe 
a sudden transition from small to large 02 values. In this 
narrow transition region, we expect two types of events: 
One with a large 02 (one big projectile residue with some 
very light fragments), the other with a small 02 (two or 
more similarly sized fragments). Events with intermedi- 
ate values of 02 should be rare. 

In order to verify whether bimodality is a ' smoking 
gun ' signal for a first order phase transition in a finite 
system, we have performed numerical simulations with 
one of the dynamical models which has frequently been 
used to interpret the multifragmentation observables, the 
Quantum Molecular Dynamics (QMD) approach 0, [1]. 
This approach simulates the entire heavy ion reaction, 
from the initial separation of projectile and target up to 
the final state, composed of fragments and single nucle- 
011s. Here, nucleons interact by mutual density depen- 
dent two body interactions and by collisions. The two 
body interaction is a paramctrization of the Briickner 
G-Matrix supplemented by an effective Coulomb inter- 
action. For this work, we have used a soft equation of 
state. The initial positions and momenta of the nucleons 
are randomly chosen, and respect the measured rnis ra- 
dius of the nuclei. Colli sions take place if two nucleons 
come closer than r = y^a/ir, where a is the energy de- 
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pendent cross section for the corresponding channel (pp 
or pn). The scattering angle is chosen randomly, respect- 
ing the experimentally measured da/d^l. Collisions are 
Pauli blocked. For details we refer to ref. 0|- For the 
later discussion it is of importance that, even for a given 
impact parameter, two simulations are not identical, be- 
cause the initial positions and momenta of the nucleons 
as well as the scattering angles are randomly chosen. 
Fig. [1] shows the INDRA experimental results (left 
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FIG. 1: (Color online) Top: most probable value of 02 in the 
quasi-projectile angular range 0cm < 90° (a^^) as a function 
of e1\2j the total transverse energy E± 12 in the quasi-target 
angular range dcm > 90°, scaled by Ecm, the energy per nu- 
cleon of the system in the center of mass . We display INDRA 
exprimental results (left panels) extracted from [3| and QMD 
simulations (right panels) for the Au-|-Au collisions at three 
different bombarding energies. Bottom: differential reaction 
cross-section (linear color scale, arbitrary unit) of a^^as a 
function of E'f^-^^, in the transition region, for Au+Au at 
lOQA.AieV bombarding energy. We show the INDRA experi- 
mental data and the filtered QMD simulations (left and right 
panels respectively). As in for calculating 02, in both ex- 
perimental and QMD results, it is required that at least 80% 
of the total charge of the projectile is detected by the INDRA 
set-up in the forward hemisphere. 



panels) in comparison with those of QMD calculations 
(right panels). The calculations have been acceptance 
corrected. Qualitatively, we see the same behavior also 
in the unfiltercd simulations. In the top row, we present 
the most probable value of 02 in the quasi-projectile as a 



function of E±_ 12 /{Ec.m./A) in the quasi-target for dif- 
ferent beam energies, where Ec.m./A is the center of mass 
energy per nucleon of the colliding system. Accordingly 
with [l| , the two observables have been determined in dis- 
tinct angular ranges in order to minimize possible correla- 
tions between the total transverse kinetic energy of light 
particles and the size of the two biggest fragments in- 
side the same spectator (quasi-projectile or quasi-target). 
We observe that in the experiment as in the calculation, 
the sudden transition between large and small 02 val- 
ues scales with E^^^/ [Ec.m. / A) . Even the numerical 
value of this transition agrees between experiment and 
theory. In order to see whether this phenomenon sur- 
vives at higher incident energies, the simulations have 
been extended up to \^QA.MeV bombarding energy. We 
observe that it is the case. 

The bottom row shows the transition interval in detail. 
Here, we display the differential reaction cross-section of 
as a function of Ef\^ for Au-f Au at IQQA.MeV. 
From the experimental data, we observe that there is 
no smooth transition between the two event classes. In 
the simulations as well as in the experiment we see two 
maxima for , separated by a minimum of the dis- 
tribution. QMD simulations reproduce the experimental 
findings qualitatively and quantitatively. 

In QMD simulations, the system does not even come 
to a local thermal equilibrium. It is therefore necessary 
to explore the origin of the observed dependence of 02 as 
a function of Ei_ 12 . The first step toward this goal is 
to identify when, in the course of the reaction, the frag- 
ment pattern is determined. This is all but trivial. Frag- 
ments can easily be identified at the end of the heavy ion 
reaction, when they are clearly separated in coordinate 
space, by a minimum spanning tree procedure. At ear- 
lier times, however, they overlap in coordinate space and, 
consequently, another method has to be employed. It has 
been proposed by Dorso and Randrup , and later been 
verified in QMD simulations that an early identifi- 
cation of fragments is possible if one uses in addition the 
momentum space information: If one identifies at each 
time step during the simulation the most bound con- 
figuration, one can establish that the fragment pattern 
changes only little during the time, and that the early 
identified fragments are the prefragments of the finally 
observed fragments. The most bound configuration in a 
simulation with N fragments is that in which 
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is minimal. Ei is the binding energy of the fragment i 
which contains m(i) nucleons and is given by 
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{Pk- < Pr >Y + V Vkl 
k=l k<l 

where < > is the average momentum of the nucleons 
entrained in the fragment i. Please note that Ebmd does 
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not contain the interaction among fragments. Therefore, 
its numerical value can vary, although the total energy 
is conserved in the simulations. The most bound config- 
uration has to be determined by a simulated annealing 
procedure With this procedure, the fragment multi- 
plicity of a given event can already be determined when 
the system has passed the highest density and is starting 
to expand. Later, the prefragments may still emit some 
nucleons, but the nucleons which are entrained at the 
end in the fragment are part of this prefragment. These 
methods allow us to trace back at which time point in 
the reaction it is determined whether the event has a 
large or a small final a2 value. Fig. [2] shows QMD re- 
sults for Au+Au at 80, 100 and IbOA.MeV bombarding 
energy. We see that, whatever the incident energy, the 
two event classes are already formed shortly after the 
system has passed the highest compression stage. There 
we observe the highest rate of hard NN collisions. These 
collisions transport nucleons in unoccupied regions of mo- 
mentum space leaving behind holes in momentum space 
which create in time (due to the different trajectories) 
holes in coordinate space. These holes weaken the bind- 
ing of projectile and target matter, leading to a fragment 
formation still at rather high density. Within statisti- 
cal models [1, , droplets are created with a spherical 
shape and have a normal density. This is only possible at 
a density of the system of p < po/3. The observation of 
an early creation of cluster partitions, above the critical 
density, as been found too in jTl| with lattice-gas calcula- 
tions where dropplcts are also defined according to energy 
considerations. The same conclusion has been obtained 
in from classical molecular dynamics (CMD) calcu- 
lations where particles interact through Lennard-Jones 
plus Coulomb potentials. 

Consequently, bimodality in QMD has nothing to 
do with the final state interaction, or with how the 
neck between projectile and target residue finally breaks. 
Whether we find a multifragmentation or a heavy residue 
event is determined when projectile and target nucleons 
still overlap almost completely in coordinate space Q. 
One may conjecture that, due to the random character 
of the scattering angle, events with the same E± 12 de- 
celerate differently, and, therefore, a different behavior of 
the average momenta may be at the origin of the differ- 
ent a2 values. For this purpose, we study with Au+Au 
at 15QA.MeV incident energy, at 60 fm/c, when 02 is 
decided, the dependence with the final a2 of the average 
momentum of all target nucleons which are at the end 
entrained in A > 4 fragments. In fig. [3l we display their 
longitudinal and transverse momentum as a function of 
a2. Clearly, both average momentum are almost inde- 
pendent of the final 02. The fiuctuations of the momenta 
around the mean values are by far larger than the dif- 
ference between the mean values for different 02 values. 
This excludes mean deceleration of the simulation events 
as reason for the different reaction scenarios. 

Obviously fiuctuations around the mean values are at 
the origin of the different event classes. This phenomenon 
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FIG. 2: (Color online) Results of QMD simulations of Au+Au 
at 80, 100 and 150 A.MeV incident energies (left, middle and 
right panels respectively). Top: function of time for 

the two different event classes: final (at 200 fm/c) a^^< 0.4 
corresponding to multifragmentation events and final a^^> 
0.4 corresponding to a projectile residue. Bottom: central 
density of the system as a function of time. 



is known from nonlinear theory [14| , and is called "bifur- 
cation" . We see here in a system with a very limited par- 
ticle number, a nonlinear behavior. Can we understand 
where it comes from? In order to answer this question, 
we go back to the search of the most bound configura- 
tion. For E±_ 12 values below the transition (and hence 
large 02 values), the most bound configuration is a large 
residue. Above the transition (where 02 is small), several 
small fragments give a more bound configuration. The 
sum of the internal kinetic energies of the clusters 

^ N rn(i) 
Emult = 7;—y2y2 (Pk- <Pi>f 

2m ^-^ ^-^ 

i=l fe=l 

is there smaller than 

^ m{i)+va{2) + ..+m{N) 
Eres = 7r- y] (Pk- < P >)^ 

fe=i 

the internal kinetic energy for a residue configuration, 
and compensates the increase of the attractive potential 
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FIG. 3: (Color online) QMD simulations of Au+Au at 
IbOA.MeV incident energy, at 60 fm/c: differential cross- 
section (colored contour levels, linear scaled) of the longitu- 
dinal (left) and transverse momentum (right) of all nucleons 
which are finally entrained in a fragment of size ^ > 4 as a 
function of final (at 200 fm/c) a^^. The symbols represent 
the mean values of momentum. 



energy 

^ m(l)+m(2) + ..+m(7V) N m{i) 

Vres — Vmult = ^ ( ^ ^kl ~ ^ ^ Vkl)- 

k,l=l i=l k,l=l 

In the transition region, we see that Emuit + Vmuit ~ 
Eres + Vres- In somc events, both configurations differ 
by IQQkeV only. Therefore it may happen that the scat- 
tering angle of one single nucleon-nucleon collision, which 
is - see above - randomly chosen, determines the type of 
the most bound configuration. 

It is interesting to sec the differences and similarities 
of the origin of bifurcation in a statistical model as com- 
pared to the analysis of QMD events. In both cases, the 
energy is the essential quantity. In the statistical model, 
there is, for a given number of nucleons in a given vol- 
ume, a small range of total energies for which the number 
of micro-states with one residue is of the same order of 
magnitude as the number of micro-states with many frag- 



ments. (In order to count the micro-states, it is assumed 
that the fragments are in one if their eigenstates, some- 
times parameterized by a level density formula.) In this 
energy range, bimodality appears as a global property 
of the systems which is dependent on the total energy 
of all nucleons present in the reaction. In QMD events, 
the essential quantity is the total binding energy of the 
nucleons bound in medium size or large clusters. As ex- 
plained above, in the transition region, this energy is al- 
most identical for a multifragment and for a residue con- 
figuration. Therefore, both configurations appear, and 
we see bimodality. The fragments are not in the ground 
state, their nucleons are not isotropic neither in coordi- 
nate space nor in momentum space. Thus, bimodality 
is a local quantity in QMD simulations, depending only 
on the total binding energy of a subset of the nucleons. 
Therefore, in QMD, bimodality makes no reference to a 
statistical or thermal equilibrium, neither of the system 
nor of the population of the excited states of the frag- 
ments. 

In summary, we have shown that the experimentally 
observed bimodality, the sudden transition between a 
residue and a multifragment exit states, and the existence 
of a small interval in E±_ 12 in which both channels are 
coexistent, is in quantitative agreement with the result of 
QMD simulations. Even the scaling of this transition re- 
gion with the center of mass energy of the system is well 
reproduced. From a detailed investigation of the reaction 
mechanism in QMD, we have seen that bimodality has 
properties observed in nonlinear systems: The system 
shows bifurcation as a function of the control parame- 
ter E± 12 . Fluctuations around the mean value in the 
longitudinal and transverse momentum decide which exit 
channel the simulation will take. 

Being reproduced in statistical as well as in dynamical 
models, bimodality reflects the same ambiguity already 
observed for other observables • 
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